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Abstract 

In the analysis of large random wireless networks, the underlying node distribution is almost ubiquitously 
assumed to be the homogeneous Poisson point process. In this paper, the node locations are assumed to 
form a Poisson clustered process on the plane. We derive the distributional properties of the interference 
and provide upper and lower bounds for its CCDF. We consider the probability of successful transmission 
in an interference limited channel when fading is modeled as Rayleigh. We provide a numerically integrable 
expression for the outage probability and closed-form upper and lower bounds. We show that when the 
transmitter-receiver distance is large, the success probability is greater than that of a Poisson arrangement. 
These results characterize the performance of the system under geographical or MAC-induced clustering. 
We obtain the maximum intensity of transmitting nodes for a given outage constraint, i.e., the transmission 
capacity (of this spatial arrangement) and show that it is equal to that of a Poisson arrangement of nodes. 
For the analysis, techniques from stochastic geometry are used, in particular the probability generating 
functional of Poisson cluster processes, the Palm characterization of Poisson cluster processes and the 
Campbell-Mecke theorem. 

I. Introduction 

A common and analytically convenient assumption for the node distribution in large wireless networks is 
the homogeneous (or stationary) Poisson point process (PPP) of intensity A, where the number of nodes 
in a certain area of size A is Poisson with parameter XA, and the numbers of nodes in two disjoint areas 
are independent random variables. For sensor networks, this assumption is usually justified by claiming that 
sensor nodes may be dropped from aircraft in large numbers; for mobile ad hoc networks, it may be argued 
that terminals move independently from each other. While this may be the case for certain networks, it is 
much more likely that the node distribution is not "completely spatially random" (CSR), i.e., that nodes 

Part of the material in this paper has been presented at the 2006 Asilomar conference. 
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are either clustered or more regularly distributed. Moreover, even if the complete set of nodes constitutes 
a PPP, the subset of active nodes {e.g., transmitters in a given time-slot or sentries in a sensor network), 
may not be homogeneously Poisson. Certainly, it is preferable that simultaneous transmitters in an ad hoc 
network or sentries in a sensor network form more regular processes to maximize spatial reuse or coverage 
respectively. On the other hand, many protocols have been suggested that are based on clustered processes. 
This motivates the need to extend the rich set of results available for PPPs to other node distributions. The 
clustering of nodes may be due to geographical factors, for example communicating nodes inside a building 
or groups of nodes moving in a coordinated fashion. The clustering may also be "artificially" induced by 
MAC protocols. We denote the former as geographical clustering and the latter as logical clustering. 

A. Related Work 

There exists a significant body of literature for networks with Poisson distributed nodes. In [1] the 
characteristic function of the interference was obtained when there is no fading and the nodes are Poisson 
distributed. They also provide the probabihty distribution function of the interference as an infinite series. 
Mathar et al., in [2], analyze the interference when the interference contribution by a transmitter located at 
X, to a receiver located at the origin is exponentially distributed with parameter ||a;|p. Using this model they 
derive the density function of the interference when the nodes are arranged as a one dimensional lattice. 
Also the Laplace transform of the interference is obtained when the nodes are Poisson distributed. 

It is known that the interference in a planar network of nodes can be modeled as a shot noise process. Let 
{xj} be a point process in R. Let be a sequence of independent and identically distributed random 

functions on R'^, independent of {xj}. Then a generalized shot noise process can be defined as [3] 

Y{x)=Y^l3j{x-xj) 
j 

If is the path loss model with fading, Y{x) is the interference at location x if all nodes xj are transmitting. 
The shot noise process is a very well studied process for noise modeling. It was first introduced by Schottky 
in the study of fiuctuations in the anode current of a thermionic diode and it was studied in detail by Rice 
[4], [5]. Daley in 1971 defined multi-dimensional shot noise and examined its existence when the points {xj} 
are Poisson distributed in R**. The existence of generalized shot-noise process, for any point process was 
studied by Westcott in [3]. Westcott also provides the Laplace transform of the shot-noise when the points 
{xj} are distributed as a Poisson cluster process. Normal convergence of the multidimensional shot-noise 
process is shown by Heinrich and Schmidt [6]. They also show that when the points {xj} form a Poisson 
point process of intensity A, the rate of convergence to a normal distribution is \/A. 

In [7] , How and Hatzinakos model the interference as a shot noise process and show that the interference 
is a symmetric a-stable process [8] when the nodes are Poisson distributed on the plane. They also show 
that channel randomness affects the dispersion of the distribution, while the path-loss exponent affects the 
exponent of the process. The throughput and outage in the presence of interference are analyzed in [9]-[ll]. 
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In [9] , the shot-noise process is analyzed using stochastic geometry when the nodes are distributed as Poisson 
and the fading is Rayleigh. In [12] upper and lower bounds are obtained under general fading and Poisson 
arrangement of nodes. 

Even in the case of the PPP, the interference distribution is not known for all fading distributions and all 
channel attenuation models. Only the characteristic function or the Laplace transform of the interference 
can be obtained in most of the cases. The Laplace transform can be used to evaluate the outage probabilities 
under Rayleigh fading characteristics [9], [13]. In the analysis of outage probability, the conditional Laplace 
transform is required, i.e., the Laplace transform given that there is a point of the process located at the 
origin. For the PPP, the conditional Laplace transform is equal to the unconditional Laplace transform. To 
the best of our knowledge, we are not aware of any literature pertaining to the interference characterization 
in a clustered network. 

[14] introduces the notion of transmission capacity, which is a measure of the area spectral efficiency of 
the successful transmissions resulting from the optimal contention density as a function of the link distance. 
Transmission capacity is defined as the product of the maximum density of successful transmissions and 
their data rate, given an outage constraint. Weber et al., provide bounds for the transmission capacity under 
different models of fading, when the node location are Poisson distributed. 

B. Main contributions and organization of the paper 

In this work, we model the transmitters as a Poisson cluster process. To circumvent technical difficulties 
we assume that the receivers are not a part of this clustered process. We then focus on a specific transmit- 
receive pair at a distance R apart, see Fig[TJ We evaluate the Laplace transform of the interference on the 
plane conditioned on the event that there is a transmitter located at the origin. Upper and lower bounds 
are obtained for the CCDF of the interference. From these bounds, it is observed that the interference is 
a heavy-tailed distribution with exponent 2/a when the path loss function is ||a;||^". When the path-loss 
function has no singularity at the origin (i.e., remains bounded), the distribution of interference depends 
heavily on the fading distribution. Using the Laplace transform, the probability of successful transmission 
between a transmitter and receiver in an interference-limited Rayleigh channel is obtained. We provide a 
numerically integrable expression for the outage probability and closed-form upper and lower bounds. The 
clustering gain G{R) is defined as the ratio of success probabilities of the clustered process and the PPP 
with the same intensity. It is observed that when the transmitter-receiver distance R is large, the clustering 
gain G{R) is greater than unity and becomes infinity as R ^ oo. The gain G{R) at small R depends on 
the path loss model and the total intensity of transmissions. We provide conditions on the total intensity 
of transmitters under which the gain is greater than unity for small R. This is useful to determine when 
logical clustering performs better than uniform deployment of nodes. We also obtain the maximum intensity 
of transmitting nodes for a given outage constraint, i.e., the transmission capacity [12], [14], [15] of this 
spatial arrangement and show that it is equal to that of a Poisson arrangement of nodes. We observe that in 
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Fig. 1. Illustration of transmitters and receivers. Cluster density is 1. Transmitter density in each cluster is 3. Spread of each 
cluster is Gaussian with standard deviation a = 0.25. Observe that the intended receiver for the transmitter at the origin is not 
a part of the cluster process. The transmitter at the origin is a part of the cluster located around the origin. 

a spread-spectrum system, clustering is beneficial for long range transmissions, and we compare DS-CDMA 
and FH-CDMA. 

The paper is organized as follows: in Section [Til we present the system model and assumptions, introduce 
the Neyman-Scott cluster process and derive its conditional generating functional. In Section IlIII we derive 
the properties of interference, outage probability and the gain function G{R). In Section llVl we derive the 
transmission capacity of the clustered network. 

II. System Model and Assumptions 

In this section we introduce the system model and derive some required results for the Poisson cluster 
process. 

A. System model and notation 

The location of transmitting nodes is modeled as a stationary and isotropic Poisson cluster process (fi 
on M^. The receiver is not considered a part of the process. See Figured] Each transmitter is assumed to 
transmit at unit power. The power received by a receiver located at z due to a transmitter at x is modeled 
as hxg{x — z), where hx is the power fading coefficient (square of the amplitude fading coefficient) associated 
with the channel between the nodes x and z. We also assume that all the fading coefficients are independent 
and are drawn from the same distribution. We will sometimes use h to denote a random variable that is 
i.i.d with the power fading coefficients. Let {o} denote the origin (0,0). We assume that the path loss model 
g{x) : \ {o} ^ M+ satisfies the following conditions. 
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1) g{x) is a continuous, positive, non-increasing function of ||a::|| and 

g{x)dx < 00, Ve > 



/ 



lM?\B(o,e) 

where B{o, e) denotes a ball of radius e around the origin. 

2) 

lim ^^=1, Vj/eR2 (1) 

||a:|Hoc g(x - y) 

g{x) is usually taken to be a power law in the form ||a;||^", (1 + ||a;||")^^ or min{l, To satisfy 

condition 1, we require a > 2. The interference at node z on the plane is given by 

hi^) = ^h,9ix-z) (2) 

The conditions required for the existence of /<^(-z) are discussed in [3]. Let W denote the additive Gaussian 
noise the receiver. We say that the communication from a transmitter at the origin to a receiver situated at 
z is successful if and only if 

W + I^\[,y{z) 

or equivalently, 



> T (3) 



hgiz) ^ T 



W + l4,{z) - 1 + T 

For the calculation of outage probability and transmission capacity, the amplitude fading \/h^ is assumed 
to be Raylcigh with mean /i, but some results arc presented for the more general case of Nakagami-m 
fading. Hence the powers arc exponentially and gamma distributed respectively. We will be evaluating 
the performance of spread-spectrum in some sections of the paper. Even though we evaluate spread-spectrum 
systems (specifically DS-CDMA and FH-CDMA) we will not be using any power control, the reason being 
that there is no central base station. 

Notation: If lima;^oo f{x)/g{x) = C, we shall use f{x) ~ g{x) if C = 1, f{x) < g{x) if < C < 1 and 
f{x)>g{x) if l<C<oo. 

B. Neyman-Scott cluster processes 

Neyman-Scott cluster processes [16] are Poisson cluster processes that result from homogeneous indepen- 
dent clustering applied to a stationary Poisson process, where the parent points form a stationary Poisson 
process (pp = {xi,X2, ■ ■ ■} of intensity Xp. The clusters are of the form N^' = Ni + Xi for each Xi e 0p. The 
Ni are a family of identical and independently distributed finite point sets with distribution independent of 
the parent process. The complete process <f) is given by 

<t> = \J N\ (4) 

Note that the parent points themselves are not included. The daughter points of the representative cluster 
Nq are scattered independently and with identical distribution F{A) = f{x)dx, A c M^, around the 
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Fig. 2. (Left) Thomas cluster process with parameters Ap = l,c = 5 and a = 0.2. The crosses indicate the parent points. 
(Right) PPP with the same intensity A = 5 for comparison. 



origin. We also assume that the scattering density of the daughter process f{x) is isotropic. This makes the 
process (j) isotropic. The intensity of the cluster process is A = ApC, where c is the average number of points 
in representative cluster. 

We further focus on more specific models for the representative cluster, namely Matern cluster processes 
and Thomas cluster processes. In these processes the number of points in the representative cluster is Poisson 
distributed with mean c. For the Matern cluster process each point is uniformly distributed in a ball of radius 
a around the origin. So the density function f{x) is given by 

f^, ||a;||<a 

m = (5) 

I otherwise. 

In the Thomas cluster process each point is scattered using a symmetric normal distribution with variance 
around the origin. So the density function f{x) is given by 



A Thomas cluster process is illustrated in FiglT] Newman-Scott cluster processes are also a Cox processes [16] 
when the number of points in the daughter cluster are Poisson distributed. The density of the driving random 
measure in this case is 

■^(y) ^ f{y - x) 

Let i?o(.) denote the expectation with respect to the reduced Palm measure [16], [17]. It is basically the 
conditional expectation for point processes, given the there is a point of the process at the origin but without 
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including the point. Let v{x) : M? [0, 1] and J^j |1 — v{x)\dx < oo. When ^ is Poisson of intensity A, the 
conditional generating functional is 

= exp ^-A [1 - v{x)]dx^ (6) 
The generating functional G{v) = E (^Hxecji ^(•'^)) °f Neyman-Scott cluster process is given by [16], [18] 



G{v) = exp f-Ap f \i-m{ [ v{x + y)f{y)dy 



dx 



where M{z) = J^^oPn^"^ is the moment generating function of the number of points in the representative 
cluster. When the number of points in the representative cluster is Poisson with mean c, as in the case of 
Matern and Thomas cluster processes, 

M{z) = exp(-c(l - z)). 
The generating functional for the representative cluster Gc{v) is given by [18], [19] 

Gc{v) = M (^^^ v{x)f{x)dx^ 
The reduced Palm distribution Pq of a Neyman-Scott cluster process (j) is given by [16]-[18], [20] 

Po=P*^o (7) 

where P is the distribution of (f>, and Q'q is the reduced Palm distribution of the finite representative cluster 
process Nq. denotes the convolution of distributions, which corresponds to the superposition of (p and 
Nq. The reduced Palm distribution Q'q is given by 

^o{Y) = -_e[ J2 lr(</'-x\{0})) (8) 

xeNo 

where (j)^ = 4> + x, is a translated point process. We require the following lemma to cvahiate the conditional 
Laplace transform of the interference. Let Q{v) denote the conditional generating functional of the Neyman- 
Scott cluster process, i.e., 

We will use a dot to indicate the variable which the functional is acting on. For example Q{v{- — y)) = 

Lemma 1: Let < v{x) < 1. The conditional generating functional of Thomas and Matern clustered 
processes is 

g{v) = G{v) [ Gciv{- - y))f{y)dy. 
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Proof: Let =Y + x. From ([5]), we have 

^o(^) = E (10) 

Let f2() denote the probabihty distribution of the representative cluster. Using the Campbell-Mecke theo- 
rem [16], we get 

f lr^{4,)nUdq^)f{x)dx (11) 

Here M denotes the space of locahy finite and simple point sequences [16] on K^. Since the representative 
cluster has a Poisson distribution of points, by Slivnyak's theorem [16] we have i^'j-i-) = ^i-)- Hence 

^oiY) ^11 lyJ0)r!(d0)/(x)dx 



niY^)fix)dx (12) 
For notational convenience let ip denote Nq. Let ipy = i/j + y. Using ([7]), we have 

Giv) = [ [ U vix)Pid<P)Cll{di>) 

l[v{x)P{dc^) I n^(^)f^oW) 

G{v) [ n v{x)n[,{d^) (13) 

Giv) [ Yl v{x) [ nidijy)fiy)dy 

G{v) [ [ Yl vixMd^y)fiy)dy 



(a) 



G{v) f f l[vix~yMdi>)fiy)dy 
G{v) [ GM—V))imv 



ib) 



(a) follows from p^ . and (b) follows from the definition of G{ )- 
So from the above lemma, we have 



g{v) = expf-Ap/" [l-A/f/ v{x + y)f{y)dy) 

x[ m( [ vix-y)fix)dx)f{y)dy (14) 



The above equation holds when all the integrals are finite. Since jix) = f{—x), then vix + y)f{y)dy = 
/r2 v{x - y)f{y)dy ^v* f,so 

g{v) = expf-Ap/ \l-Miiv*f)ix))]dx) f M{(v * f){y))f{y)dy (15) 
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Likelihood and nearest neighbor functions of the Poisson cluster process, which involve similar calculations 
with Palm distributions are provided in [21]. One can obtain the nearest-neighbor distribution function of 
Thomas or Matern cluster process as D{r) = ^(ls(o,r)<^ (•))• some cases the number of points per cluster 
may be fixed rather than Poisson. The conditional generating functional, for this case is given in Appendix 

m 

III. Interference and Outage Probability of Poisson Cluster Processes 

In this section, we first derive the characteristics of interference in a Poisson clustered process conditioned 
on the existence of a transmitting node at the origin. We then evaluate the outage probability for a transmit- 
reeeive pair when the transmitters are distributed as a Neyman-Scott cluster process, with the number of 
points in each cluster is Poisson with mean c and density function f{x). 



A. Properties of the Interference I(j,{z) 

Let Ch{s) denote the Laplace transform of the fading random variable h. 
Lemma 2: The conditional Laplace transform of the interference is given by 



Proof: From ([2]) we have 



i^i4,{z){s) ^ Q {Ch{sg{- - z))) 

x£<j> 

eyi^{-shxg{x - z)) 



(16) 



Yl Ch{sg{x - z)) 



(17) 



where (a) follows from the independence of and follows from Q. ■ 
We observe from Lemma [2] and (fT5|) . that the conditional Laplace transform of the interference £/_^(z)(s) 
depends on the position z. This implies that the distribution of the interference depends on the location 
z at which we observe the interference. This is in contrast to the fact that the interference distribution is 
independent of the location z when the transmitters are Poisson distributed on the plane [9], [12]. This is due 
to the non-stationarity of the reduced Palm measure of the Neyman-Scott cluster processes. If one interprets 
Iiji{z) as a stochastic process, it is then a non stationary process due to the above reason. 

Let JCn{B) denote the reduced n-th factorial moment measure [16], [18] of a point process -0, and let 
B = Bi X . . . X Bn-i, B, e ]R2. 



ICn{B) 



X! ^b{xi,. 

Xl,...,Xn-l€ll' 



1 Xn—1 ) 



(18) 
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K^2{B{0, R)), for example, denotes the average number of points inside a ball of radius R centered around the 
origin, given that a point exists at the origin. First and second moments of the interference can be determined 
using the second and third order reduced factorial moments. The average interference (conditioned on the 
event that there is a point of the process at the origin) is given by 

= E[h]X [ g{x ^ z)IC2{dx) 
Since the process (p is stationary, JC2{B) can be expressed as [16], [22] 

where p'^'^\x) is the second order product densitjj^. So we have 



(19) 



A 

Example: Thomas Cluster Process. In this case, from [16] 

p(2)(^) 1 



E^M^)] = ^ / g{x-z)p^^\x)dx (20) 



A2 



i (~\\x\\ \ 



where A = ApC. We obtain 



E[[mz)] = £;/p,i + ^ g{x - z) exp dx (21) 

Where EI,^, is the average interference seen by a receiver located at when the nodes are distributed 
as a PPP with intensity A. The above expression also shows that the mean interference! is indeed larger 
than for the PPP. One can also get the above from the conditional Laplace transform in Lemma [2] and 
using ^'^[/^(z)] = —-j^Ci^(^z)is)\s=o- In the following theorem we provide bounds to the tail probability of 
the interference /^(z) for any stationary distribution (j) of transmitters. We adapt the technique presented 
in [15] to derive the tail bounds of the interference. We denote the tail probability (CCDF) of Icj,{z) by 
Fi{y)=¥{I^{z)>y). 

Theorem 1: When the transmitters are distributed as a stationary and isotropic point process 4> of intensity 
A with conditional generating functional Q and second order product density p^^-* , the tail probability Fi (y) 
of the interference at location z, conditioned on a transmitter present at the origij^ is lower bounded by 

^Intuitively, this indicates the probability that there are two points separated by For PPP, it is p'^' (x) = independent of 
X. Also the second order product density is a function of two arguments i.e., p^^\xi^, X2). But when the process is stationary, 
p''^' depends only on the difference of its arguments i.e., p^^^ {xi,X2) = iy{xi — X2) for all xi,X2 S M^. Furthermore if (j> is 
motion-invariant, i.e., stationary and isotropic, then v depends only on \\xi — X2\\ [16, pg 112]. 

^Note that for g{x) = \\x\\~°', E'q[I^{z)] is diverging. 

''We do not include the contribution of the transmitter at the origin in the interference. This is because the transmitter at 
the origin is the intended transmitter which we focus on. 
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Fj{y) and upper bounded by Ff{y), where 

^JC' = to)) 

Ff{,,) = l-(l-v{y))g{Fui^jJ—^'^ (23) 

where (x) denotes the CDF of the power fading coefficient h and 

2^ f ry/a{x-z) 
V{y) = ^ 9{x - z)p^''\x) vdFh{v)dx. 

Proof: The basic idea is to divide the transmitter set into two subsets (py and t/i^ where, 

(t)y = {x e 4>, Kgix - z)>y} (24) 
(jfy = {a; e (j), Kg{x -z)<y} (25) 

4>y consists of those transmitters, whose contribution to the interference exceeds y. We have I^{z) = I^^ {z) + 
I^a[z), where I(j,y{z) corresponds to the interference due to the transmitter set and l4,^{z) corresponds 
to the interference due to the transmitter set ^y. Hence we have 

Fi{y) = ni<i>y{z) + i4>i{z)>y) 

> ¥{I^^{z)>y) 
= l-P(/^^(z)<y) 

= 1 - F{(l)y = 0). (26) 

We can evaluate the probability F{^y = 0) that ^y is empty using the conditional Laplace functional as 
follows: 



-^0 n {^h^g{x-z)<y) 

= Ell[F^ 



V 



where (a) follows from the independence of h^- To obtain the upper bound 



g{x - z) 

) ) ' (27) 



Fi{y) = nh >y\hy> y)F\{y) + >y\hy< - F\{y)) 
i-q(f. 



9{--z)JJ ..-^^-^y ■'\9{--z) 

1 - (1 - P(7^ > 2/ 1 < y))g (^F^ (^(7^)) (28) 
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where (a) follows from the lower bound we have established. To evaluate P(/0 > y \ I^,^ < y) we use the 
Markov inequality (the Chebeshev inequality can also be used but is more difficult to be evaluated in this 
particular setting). We have 

>y\hy<y) = nh >y\^y=^) 
(«) 4 I 0, ^ 0) 



2 ry/y{x-z) 

y ^^j. Jo 



^h^g{x-z)<y 
f-y/aix-z) 



•■y/g{x-z) 



(h\ 1 r ry/9{x~z) 

= g{x^z) vAFh{v)p(^\x)dx (29) 

(a) follows from the Markov inequality, and (6) follows from a procedure similar to the calculation of the 
mean interference in (j20p . ■ 
In the proof of Lemma[31 we show ip{y) ~ ^22/~^^" when g{x) = This indicates the tightness of the 

bounds for large y. Lemma [3] shows that the interference is a heavy-tailed distribution with parameter 2/ a 
when the nodes are distributed as a Neyman-Scott cluster process. 

Lemma 3: For g{x) — the lower and upper bounds to CCDF Fi{y) of the interference at location 

z, when the nodes are distributed as a Neyman-Scott cluster process scale as follows for y — > oo. 

F\[y) ^ e,y-''" (30) 

F^y) ~ (^^1+^2)?/-'/" (31) 

where Oi = 7rc[(/ * f){z) + Ap] v'^/°'dFh{v) and 62 = 26*1/(0! - 2). 

Proof: See Appendix [X] ■ 
Remarks: 

1) Observe that 6*1 = ttA-V^^H^) /o°° J^^/"dFft(i/). A similar kind of scaling law with 6*1 = TT\-lp^^^\z)Eh[v^l°'\ 
and 62 = 2^1 /(a — 2) can be obtained when the transmitters are scattered as any "nice'[^ stationary, 
isotropic point process with intensity A and second order product density /^'^^(x) 7^ at x = z. 

2) A similar heavy-tailed distribution with parameter 2/ a was obtained for Poisson interference in [1], 
[15]. Since 2/ a < 1, the mean and hence the variance diverge. This can also be inferred from pT|) and is 
due to the singularity of the channel function g{x) = at the origin. For Matern cluster processes 

^ We require the conditional generating functional to have a series expansion with respect to reduced n-th factorial moment 
measures of the reduced Palm distribution [22] similar to that of the expansion of generating functional [16, p. 116] and [23]. The 
proof of the existence and the series expansion of the conditional generating functional with respect to reduced n-th factorial 
moment measures, would be of more technical nature following a technique used in [23]. If such an expansion exists it is 
straightforward to prove the scaling laws for the CCDF of interference similar to Lemma[3l with 61 = ttX~^ p^'^\z)Efi[u'^/°'] and 
6*2 = 26»i/(a - 2). 



DRAFT 



February 1, 2008 



13 



(/ * f){^) — Oi for ||z|| > 2a and for Thomas cluster processes (/ * f){z) is a Gaussian with variance 
2(7^. Hence for large z, we observe that the constants 6i become similar to that of the unconditional 
interference. This is because, the contribution of the cluster at origin becomes small as we move far 
from the origin. 

3) When the path loss function is g{x) = (1 + ||a;||")^^, the distribution of the interference more strongly 
depends on the fading model. Using a similar proof as in Lemma [3l one can deduce an exponential 
tail decay when g{x) = (1 + and Rayleigh fading. Similarly if the power fading coefficient 



o=0.25,A, =2,c=3,R=0.3,a=4 
p 




100 200 300 400 

y 



500 



600 



different fading. They were generated using Monte-Carlo simulation. Curves #1 and #2 correspond to g{x) = (1 + 
Curve #1 corresponds to Rayleigh fading and exhibits an exponential decay. Curve #2 for which h is distributed as generalized 
Pareto with parameters k = 1,6 = 0,(Tp = I (a, hypothetical power fading distribution which exhibits power law decay) exhibits 
a power law decay. Curves #3 (generalized Pareto ) and #4 (Rayleigh) correspond to g{x) = and exhibit a heavy tail 

for both fading distributions. 



follows a power-law distribution with exponent fc, the tail of the interference shows a power-law decay. 
This is because of the presence of the term y^^/a ^^"[1 — Ffi{u)]{u — y)^/"^^du in the proof. So when 
using non-singular channel models, the interference has a more intricate dependence on the fading 
characteristics rather than a simple dependence on i^/Jj/^/"] as in the singular case. This behavior 
is well understood for Poisson and unconditional Poisson cluster shot noise process [24], [25]. The 
properties of interference for different path loss models with no fading, when the nodes are uniformly 
distributed are discussed in [26]. 
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B. Success probability: ¥{success) 

Let the desired transmitter be located at the origin and the receiver at location z at distance R — \\z\\ 
from the transmitter. With a slight abuse of notation we shall be using R to denote the point (i?, 0). The 
probability of success for this pair is given by 

hg{z) 



We now assume Rayleigh fading, i.e., the received power is exponentially distributed with mean /i. So we 
have 

P(success) = / e-f^'^^sizUFiW + I^\{o}{z) < s) 
Jo 

= ^i^iz){tJ'T/g{z))£w{lJ'T/g{z)) , (33) 

When hr^ is Rayleigh we have 

U{sg{x-z)) = ——^ ^ (34) 

fj, + sg[x — zj 

At s — iiT/g(R) we observe that the above expression will be independent of the mean of the exponential 
distribution /i. 

Lemma 4- [Success probability] The probability of successful transmission between the transmitter at the 
origin and the receiver located at z g R^, when W = (no noise), is given by 



success 



)=exp|-Ap/ 1 - exp(-c/3(z,?/)) dyj 



X / eM-cPiz,y))fiy)dy (35) 



T2 

where 

Jm? ^ + g{x - y ~ z) 

Proof: Follows from ([34| and Lemma [21 ■ 
The success probability, when the number of nodes in each cluster is fixed is given in the Appendix [B) See 
Figure |4] for comparison. When the fading is Nakagami-m, the probability of success is evaluated in the 
Appendix ICl for integer m. 
Remarks: 

1) The term Ti in ()35|) captures the interference without the cluster at the origin (i.e., without condi- 
tioning); it is independeniQ of the position z since the original cluster process is stationary (can be 
verified by change of variables yi = y + z). The second term T2 is the contribution of the transmitter's 

^By this we mean the unconditional interference distribution which leads to this term does not depend on the location z. The 
term Ti does depend on g{z). 
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Fig. 4. Comparison of P(success) when the number of points in a cluster are fixed and Poisson distributed with parameter c. 

cluster; it is identical for all z with ||z|| = R since / and g are isotropic. So the success probability 
itself is the same for all z at distance R. This is because the Palm distribution is always isotropic when 
the original distribution is motion-invariant [16]. Hence we shall use (3{R,y) to denote P{z,y) where 
z — i?e'^. We shall also use R and {R, 0) interchangeably and will be clear by the context. 

2) From the above argument we observe that P(success) depends only on ||z|| = i? and not on the angle 
of z. So the success probability should be interpreted as an average over the circle ||z|| = R, i.e., the 
receiver may be uniformly located anywhere on the circle of radius R around the origin. For large 
distances R, there is a very high probability that the receiver is located in an empty space and not in 
any cluster. Hence for large R the success probability is higher than that of a PPP of the same intensity. 
If the receiver is also conditioned to be in a cluster, we have to multiply (at least heuristically) by a 
term that is similar to T2 and this would significantly reduce the success probability. 

3) From LemmalU we have P(success) — i?5[exp(— s/0(z))] evaluated at s = /iT/g(z). If fiT/g{z) is small, 
and J g{z)dz < 00 (i.e, finite average interference) then, P(success) < Pp{X). This follows from (pij) 
and the fact that E'q[Ici,{z)] is the slope of the curve E'Q[exp{—sI^{z))] at s = 0. This implies that 
at small distances, spread spectrum (DS-CDMA) works better with a Poisson distribution of nodes. 
(If the distance R is large, then the spreading gain has to increase approximately like g{z) to keep 

/ g{z) small.) 

4) Let M be the DS-CDMA spreading factor. We have P(outage) = V{I^{z) > ^ha:g{z)). For g{x) = 
||a;||^", we have the following scaling law for the outage probability with respect to the spreading gain. 

(a) (b) 

l9lJ?2M'2/ay2/a^^j^-2/al < p^Q^jt^ge) < —±-eiRH'I-'^^°'T^/'^Eh[i'~'^^°'], (37) 

a — 2 
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where Eh[i^^^^"] — i'^^^°'dFi^2(iy). (a) and (b) follow from Lemma [31 Also observe that these scaling 
bounds arc valid for any fading distribution for which iJ/Ji^^^/"] < oo. Similar scaling laws with the 
exponent of M being —Ija can be obtained when the transmitters are Poisson distributed on the 
plane. When the fading is Rayleigh i.e., h ^ exp(/i), the lower bound is 

7rc[(/ * + Ap]r(l + ^)r(l - ^)i?2M-2/"r2/a 

and the upper bound is q;/(q;— 2) times the lower bound. r(z) represents the standard Gamma function. 
We now derive closed form upper and lower bounds on P(succcss). 
Lemma 5: [Lower bound] 

P(success) > Pp{X)Pp{cf*) (38) 

where -Pp(A) denotes the success probability when is a PPP, /* = supy^-g2{f * f)(y), and A = ApC. 

Proof: The first factor in ([35|l , Ti can be lower bounded by the success probability in the standard PPP 
Pp(A), and the second factor can be lower bounded by Pp{cf*). From (I35p and the fact that 1 — exp(— Ja;) < 
6x, 5 > 0, we have 

P(success) > exp ( - ApC / P{R, y)dy) (39) 

Tcrml 

X / exp{-cP{R,y))f{y)dy 



Torm2 



Terml = cxp ( - A / f3{R, y)dy 
exp (-a/ 



dy 



T 

PpW (40) 



(a) follows from change of variables, interchanging integrals and using J f{x) — 1. 

Term2 = / exp(-c/3(i?, y))/(y)dy 
Since exp(— x) is convex and f{x) > 0, J f{x) — 1, Using Jensen's inequality {Ef{x) > f{E{x))) we have, 

Term2 > exp ( - c f (3{R, y)f{y)dy ) 
Changing variables and using j[x) — /(—a;), we get, 

Term2 > cxp ( ~ c [ J^""^ [ f{x + z- y)f{y)dydx ) 

> exp(-c[ ^J^''\ (f*f){x + z)dx) (41) 
^ Jr2 £i£ii + „(x) ' 



Hence 

Term2 > Pp(c/*) (42) 
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Since f £ Lp , by Young's inequality [27] we have /* < ||/||p||/||<j, where 1/p+l/q = 1 (conjugate exponents). 
For a > I/Vtt (Matern) and a > l/\/27r (Thomas), we get P(success) > Pp{X)Pp{c). In general, /* < 
||/||oo||/||i, which is I/tto^ for Matern and l/27r(T^ for Thomas processes. In the latter case, when / is 
Gaussian, / * / is also Gaussian with variance 2(7^, hence /* < l/Ana^. From [9], we get (by change of 
variables) : 



Pp(A) =expf-A [ P{R,y)dy 



(43) 



We have 

. for g{x) = \\x\\-°', Pp{X) = cxp(-Ai?2r2/"C(a)) [9], where C{a) = (27rr(2/a)r(l - 2/a))/a = 
2ficsc(2Va). 

. for gix) = (1 + 11x11")-!, Pp{X) = exp(-ArC(a)(T + 5(i?))2/"-i5(i?)-2/a). 

Let Pj = Jjgs /9(^, J/)dj/, P = supj,gjj2 /3(i?, y) and / = supygj|2 /(y). By Holders inequality we have /3 < 
min{l,//3/(i?)}. Also let k = l3{R,y)f{y)dy. 
Lemma 6: [Upper bound] 

P(success) < Pp i — ^ I (44) 
" Vl + c/3/ 

Proof: Neglecting the second term T2 and using exp(— (5a;) < 1/(1 + 5x), we have 



^(success) < exp(-A,|^Jl-^^^l^]d,) 



V Jg_2 1 + c/3(i?,y) 
< exp(--^ / l3{R,y)dy) (45) 



From the above two lemmata, we get 

Pp{X)Pp{cr) < P(success) < Pp (^^] (46) 

V 1 + c/3 / 

from which follows P(success) —> ^'p(A) as — > as expected. In Lemma IHl we have neglected the 
contribution of the transmitter's cluster. We derive the following upper bound in the proof of Lemma [H 

P(success) < Pp(A)exp (A/3/zy(c/3)) 1 - (l - i/(c/3)) ck (47) 

where ^{x) — (exp(— x) — 1 + x)/x. Substituting for v{x), we have 

'a(1 -exp(-c/3))' 



""(success) < Pp 



1 - ( 1 - exp(-c/3) ) I 



(48) 



c/3 

is a tighter bound than the bound in Lemma [6l but not easily computable due to the presence of k (for 
a given R, T and cr, k and (3* are constants). In ((48|) . the outage due to the interference by the transmitting 
cluster is also taken into account. 
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The proof of Lemmata [S] and [5] also indicates that it is only by conditioning on the event that there 
is a point at the origin that the success probability of Ncyman-Scott cluster processes can be lower than 
the Poisson process of the same intensity. This implies that the cluster around the transmitter causes the 
maximum "damage". So as the receiver moves away from the transmitter, the Neyman-Scott cluster process 
has a better success probability than the PPP. So, it is not true in general that cluster processes have 
a lower success probability than PPPs of the same intensity. For example from Figure [5l we see that for 



a=4, T.0.1, a.0.6 

0.8 I 1 1 1 r- 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

R 



Fig. 5. Comparison of success probability for cluster and Poisson process of intensity 2 

R < 0.8, the PPP has a better success probability than the Matern process. In Subsection IIII-CI we give a 
more detailed analysis, which reveals that a PPP with intensity XpC has a lower success probability than a 
clustered process of the same intensity for large transmit-receiver distances. On the other hand, for small R, 
the success probability of the PPP is higher. 

C. Clustering Gain G{R) 

In this subsection we compare the performances of a clustered network and a Poisson network of the 
same intensity with Rayleigh fading. We deduce how the clustering gain depends on the transmitter receiver 
distance. We use the following notation, 

PiiR,c,\p) = exp(-Ap/ l-eM-cl3{R,y))dy) (49) 

P2{R,c) ^ [ exp{~cPiR,y))f{y)dy (50) 

So P(success) = Pi{R,c, Xp)P2{R,c). P2 is the probability of success due to the presence of the cluster at 
the origin near the transmitter. Pi is the probability of success in the presence of other clusters. Interference 
from these other clusters contributes more to the outage when R is large. This is also intuitive, since as 
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the receiver moves away from the transmitting cluster, the interference from the other clusters starts to 

dominate. We define the clustering gain G{R) as 

^^^^ ^ P,{R,c, Xp)P2{R,c) 
Pp{XpC.) 

The fluctuation of G{R) around unity indicates the existence of a crossover point R* below which the PPP 
performs better than clustered process and vice versa. The values of G{R) at the origin and infinity indicate 
the gain of scheduling transmitters as clusters instead of being spread uniformly on the plane. So it is 
beneficial to induce logical clustering of transmitters by MAC if G{R) > 1. 

We first consider G{R) for large R, i.e., lim/j^oo G{R). By the dominated convergence theorem and ([T]), 
we have 

lim F2(i?,c) = / expf-c/" lim <ix)f{y)dy 

Tg{x-y-R) 

Also from the derivation of upper bound we have Pi{R,c, Xp) < Pp (y+e^)' Hence from the definition of 
Pp{x) we have, hmji^oo Pi{R, c, Xp) = 0. Hence for large R, 

PiiR,c,Xp) <P2iR,c) (52) 

So for large R, most of the damage is done by transmitting nodes other than the cluster in which the intended 
transmitter lies. 
Lemma 7: 

1™ u^i^-^l ^ = (53) 
Proof: See Appendix [Dl ■ 

Hence for large R, p^^j^'^g"^^ ) < exp( .^_^~^^y ). From ([5T|) we have Pp{Xpc) < Pi{R,c, Xp)P2{R,c), for large R, 

i.e., G{R) > 1 for large transmit-receive distance. We have lim/j^oo = oo. Hence the Poisson point 

process with intensity XpC, has a lower success probability than the clustered process of the same intensity for 

large transmit receiver distances. 

For small R, G{R) depends on the behavior of the path loss function, g{x) at ||a:;|l = 0. We consider the 
two cases when the channel function is singular at the origin or not. 

1) lini||j.[|^o .9(2;) = OO-' III this case we observe that G(0) = 1. But at small R, G{R) is less than 1. We 
have the following lemma. 

Lemma 8: If (/ * f){x) > \\x\\ for small \\x\\ and g(x) = |la;|l~", then for small R, 

P(success) < PpiXpc) (54) 
Proof: See Appendix [El ■ 
Note that f{x) for Matern and Thomas cluster process have the required property. Hence when g(x) = ||a;|| 
the PPP with intensity XpC, has a higher success probability than the clustered process of the same intensity 
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for small transmit receiver distance. Lemma [5] and the fact that G{oo) — oo also indicate the existence of a 
crossover point R* between the success curves of the PPP and the cluster process. So it is not true in general 
that the performance of the clustered process is better or worse than that of the Poisson process. This is 
because, for the same intensity, a clustered process will have clusters of transmitters (where interference is 
high) and also vacant areas (where there are no transmitters and interference is low), whereas in a Poisson 
process, the transmitters are uniformly spread. 

2) lini||^||_»o 5(2;) — g < 00: Pi{R,c, Xp) can be written as 



Pi{R,c,Xp) = Pp(Apc)exp 



~ 00 



^-^c"/3(i?,y)"dj/ 



>o 



Hence G(i?) can also be written as follows 

G{R) = P2{R,c)exp(^Xpcl^J{R,y)Tj{c,R,y)dy^ (55) 

where r]{c,R,y) ~ i^{c(3{R,y)), with i^{x) — (exp(— .t) — l + x)/x. Observe that < rj{c,R,y) < l,Va; > 0. If 
the total density of the transmitters is fixed i.e., X = XpC is constant, how does G{R) behave with respect 
to c? We have the following lemma which characterizes the monotonicity of G{R) with respect to c. 



Lemma 9: Given A = ApC is constant, G{R) is decreasing with c, i. 



< , Vc> iff A < X*{R,T), 



where 



X*{R,T) 



Proof: From ((55)) . 

G{R) 



2f^._f3{R,y)f{y)dy 
4./3(i?,y)2dy 



P2{R, c) exp 



XpC / l3{R,a)ri{c, R, a)da 



exp 



( - cP{R, y) + >^ I a)vic, R, a)da) f{y)dy 



(56) 



We have 

dG{R) 
dc 



|c=o = f3{R, z)/2 and '^''^ j'^'^'' is decreasing in c. 



-I3{R, y) + X I (3{R, a) ^^^i^ii^da 



dc 



exp[A / (3{R,a)ri{c, R,a)da] 



exp - c(3{R,y) + X j ^ P{R, a)ri{c, R, a)da^ f{v)dy 

exp - c(i(R, y))/(y)dy 



d?7(c, R, a) 



PiR,y) + X / PiR,z) 

iR2 dc 



T2(e) 

Since r]'{c, R, z) is decreasing in c, we have T2(c) is decreasing in c. So a necessary and sufficient condition 
for < Vc > is T2(0) < 0. We want 



dc 



^2(0) - / 



X < 



-/3(i?,2/) + - / P^{R,z)dz 

2j^,PiR,y)f{y)dy 
4. P^R,z)dz 



fiy)dy < 



(57) 



Remarks: 
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1) Since P{0,y) ^ 0, we have that, G{0) is increasing with Xp (hke exp(Ap)), and hence will be greater 
than 1 at some Xp for a fixed c. 

2) We have linig^o G{R) = 1 and specifically G(0) = 1 at c 0. 

3) From Lemma |9] and Remark 2 we can deduce G{R) < 1, Vc > if A < X*{R,T) i.e., the gain G{R) 
decreases from 1 with increasing c if the total intensity of transmitters is less than X*{R,T). 

4) Since G{R) is continuous with respect to R, G{R) is close to G(0) for small R. 

5) From Figure [71 we observe that G{R) increases monotonically with R. 

In Figure [HI A*(0,T) is plotted against T. We provide some heuristics as to when logical clustering does not 
perform better than a uniform distribution of points: 

• The exact value of R at which G{R) crosses 1 is difficult to find analytically due to the highly nonlinear 
nature of G{R). If such a crossover point exists (depends on the path-loss model) we will denote it by 
R*. 

m If g{x) = ||a;|| it is better to induce logical clustering by the MAC scheme if the link distance is larger 
than R*. Otherwise it is better to schedule the transmissions so that they are scattered uniformly on 
the plane. 

• If g{0) < oo and for a constant intensity XpC, it is always beneficial to induce clustering for long-hop 
transmissions. When R is small the answer depends on the total intensity XpC. If XpC < X* (0, T) then 
G{0) < 1 by observation 3, and hence G{R) < 1 for small R by observation 4. Also when XpC < A*(0, T), 
it is better to reduce logical clustering by decreasing c and increasing Xp, since G(0) is a decreasing 
function of c. From Figure [6]we observe that A*(0, 0.5) w 1.26 when g{x) = (1 + ||x||*)^^ and a = 0.25. 
In Figure [71 G{R) is plotted for XpC = 0.75, 9 for the same values of cr, a and the same channel function 
as of Figure [HI When XpC = 9 > A* (0,0.5), we observe that the gain curve G{R) is approximately 10 
at the origin and increases. When ApC — 0.75 < A*(0,0.5), G(R) starts around 0.25 and crosses 1 at 
i? w 1.2. We also observe that G{R), for the non-singular g{x), seem to increase monotonically. We also 
observe that the gain function for g{x) = \\x\\^" decreases from 1 initially and then increases to infinity. 

• For DS-CDMA, the value of T is smaller by a factor equal to the spreading gain. From Figure [H we 
observe that the threshold A*(0,T) for clustering to be beneficial at small distances increases with 
decreasing T. Hence for a constant intensity of transmissions XpC, the benefit of clustering decreases 
with increasing spreading gain for small link distances. So for DS-CDMA (for a large spreading gain) 
it is better to make the transmissions uniform on the plane for smaller link distances and cluster the 
transmitters for long-range communication. 

• For FH-CDMA, the total number of transmissions XpC is reduced by the spreading gain while T remains 
constant (see Figure [6]). Hence ApC < A*(0, T) for small distances and one can draw similar conclusions 
as that of DS-CDMA. The relative gain between FH-CDMA and DS-CDMA with clustering is more 
difficult to characterize analytically. 
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Fig. 6. A*(0, T) versus T for g{x) = (1 + ||a;[|^)~^, cr = 0.25. The region below the curve consists of all the pairs of (T, A = Apc) 
such that G(0) < 1. "Normal operating point" denotes a pair (T, A) that lies above the curve (T, A*(0, T)). Suppose we use FH- 
CDMA, the total intensity decreases by a factor of spreading gain and hence we move vertically downwards into the G(0) < 1 
region. If DS-CDMA is used, the threshold T decreases by a factor of spreading gain and hence we move horizontally towards 
the left into the G(0) < 1 region. 



IV. Transmission Capacity of Clustered Transmitters 

It is important to understand the performance of ad hoc wireless networks. Transmission capacity was 
introduced in [12], [14], [15] and is defined as the product of the maximum density of successful transmissions 
and their data rate, given an outage constraint. More formally, if the intensity of the contending transmitters 
is A with an outage threshold T and a bit rate b bits per second per hertz, then the transmission capacity 
at a fixed distance R is given by 

C(e, T) = 6(1 - e) sup{A : P(A, T) > 1 - e} (58) 

A 

where P(A, T) denotes the success probability of a given transmitter receiver pair. More discussion about 
the transmission capacity and its relation to other metrics like transport capacity is provided in [15]. Note 
that the results proved in [12], [14], [15] are for Poisson arrangement of transmitters. 

In this section we evaluate the transmission capacity when the transmitters are arranged as a Poisson 
clustered process. We prove that for small values of e, the transmission capacity of the clustered process 
coincides with that of the Poisson arrangement of nodes. We also show that care should be taken in defining 
transmission capacity for general distribution of nodes. For notational convenience we shall assume 6—1. For 
the clustered process, P(A, T) denotes the success probability of the cluster process with intensity A — XpC 
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G(R), a=4, a=0.25 




0.5 1 1.5 

R 



Fig. 7. G{R) versus R, a = A, a = 0.25. Observe that the gain curves #2 and #3, which correspond to the singular channel, 
start at 1 decrease and then increase above unity. For the gain curve #4, the total intensity of transmitters is 3 * 0.25 = 0.75 
which is loss than the threshold A* (0,0.5) ~ 1.26. Hence the gain curve for this case starts below unity at = and then 
increases. For the gain curve #1 the total intensity is 9 > 1.26. By chance, in the present case the gain curve #1 starts around 
10 and increases. 



and threshold T. Let Pi{X,T), Pu{\,T) denote lower and upper bounds of the success probability P{X,T) 
and the corresponding sets Ai,Au defined by := {A : P^{X,T) > 1 — e} for x G We then have 

Ai C A C Au which implies 

sup A; < sup A < SUp^ti. (59) 

Let Ci{e,T) ~ sup A; and Cu{e,T) — sup An denote lower and upper bounds to the transmission capacity. 

For a PPP we have from ([321), Pp{^iT) ~ exp(— A/3/) (/?/ does not depend on A). Hence the transmission 
capacity of a PPP denoted by Cp(e,T) is given by 



- (60) 
Pi 

For Neyman-Scott cluster processes, the intensity A = \pC. We first to try to consider both \p and c as 
optimization parameters for the transmission capacity, i.e. 

C(e, T) := (1 — e) sup{ApC : Ap > 0, c > 0, outage-constraint} (61) 

without individually constraining the parent node density or the average number of nodes per cluster. 
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Lemma 10: The transmission capacity of Poisson clustered processes is lower bounded by the transmission 
capacity of the PPP, 

C{e,T)>Ci{e,T)^Cp{e,T) (62) 
Proof: From Lemma[51 we have Pi{\, T) = Pp{Xpc)Pp{cf*). So to get a lower bound, from ([S^]) we have 
to find 

snp [XpC-.XpC + cf* < ^1^(t^)=^T^} ^^^^ 

This maximum value of XpC is attained when, Xp oo, while c ^ 0, such that cXp = Cp{e,T)(l — e)^^. So 
we have C/ (e, T) = Cp (e, T) . ■ 
Also observe that Xp oo and c — > 0. This corresponds to the scenario in which the clustered process 
degenerated to a PPP. We also have the following upper bound. 

Lemma 11: Let p{T) = k/$ with J l3{R,y)f{y)dy. For e < 1 - e-P^'^\ we have 

C{e,T) < Cu{e,T) = Cp{e,T) (64) 
Proof: See Appendix [F] ■ 
Theorem 2: For e < 1 - e~P^'^^ we have C(e,T) = Cp(e,T). 
Proof: Follows from the Lemmata [TU] and [TT] ■ 

From the above two proofs, when e is small, the transmission capacity is equal to the Poisson process of 
same intensity. This capacity is achieved when Ap oo and c — > 0. This is the scenario in which the cluster 
process becomes a PPP. This is due to the definition of the transmission capacity as C(e, T) := sup{ApC : 
Ap > 0, c > 0, outage-constraint} where we have two variables to optimize over. 

Instead we may fix Ap as constant and find the transmission capacity with respect to c. So we define 
constrained transmission capacity as 

C*(e, T) := Ap(l — e) sup{c : c > 0, outage-constraint} (65) 

We have the following bounds for C* (e, T) 
Theorem 3: 

h&2ll <C*(e,T)< , KCp{e,T) ^ ^^^^ 

Ap + r niax|0,Ap-|-ln(^j| 
Proof: From the lower bound on P(success), we have to find 



sup{c-:ApC-+c-f < ^ln(^j^^} (67) 



_ Cp{e,T) ^ 



So we have C;(e,r) = Cp(e, r)/(/* -f Ap). 

From the upper bound on P(success),we have to find 



^ l + cB I3i \l-ej 1-e ) 
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R=0.25,o=0.25,E=0.1,T=1,A =1 
p 
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a 



Fig. 8. Upper and lower bounds of C*{e,T) versus a, g{x) = \\x\\ " ,T = l,a = 0.25, e = 0.1, Ap = 1 



One can also derive an order approximation to the constrained transmission capacity when e is very smaU. 
We have the following order approximation to transmission capacity. 

Proposition 1: When Ap is fixed, the constrained transmission capacity is given by 

- <'-'>(Ai7TI + *') 

when e ^ 0. 

Proof: Let 7(c) denote the outage probability, i.e., 

7(c) = f-expf-Ap/" f-exp[-c/3(i?,y)]dj/} /" exp(-c/3(i?,2/))/(y)d2/ (70) 

We have d7(c)/dc > 0, which implies 7(c) is increasing and invertible and hence C*(e, T) = Ap(f — e)7~^(e). 
We approximate 7^^(e) for small e by the Lagrange inversion theorem. Observe that 7(c) is a smooth function 
of c and all derivatives exist. Expanding 7~^(e) around e = by the Lagrange inversion theorem and using 
7(0) = yields 

— / c " 



|e=o + o(e) 



c-=o-T (71) 



7(c) 
(a) e 



o e 



Ap/3/ + /c 

where (a) follows by applying de L'Hopital's rule. 
We have the following observations 

1) The constrained transmission capacity increases (slowly) with Ap. 
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2) We also observe that the constrained transmission capacity for the cluster process is always less than 
that of a Poisson network (see Figure [8]) and approaches Cp{e,T) as Ap — + oo. 

3) When FH-CDMA with intra-cluster frequency hopping is utilized, we have the cluster intensity c 
reduced by a factor M (spreading gain). One can easily obtain the constrained transmission capacity 
of this system to be 



C>^(e,T) = (!-£) 



€Xr,M 



When DS-CDMA is used, the constrained transmission capacity is C*jjg[e,T) — C*{e,T/M). When 
the transmitters are spread as a Poisson point process, we have from [28], [29] 

fCFH{e,Ty 



In 



= (1 - 2/a)ln(M). 



In Figure [H we plot ln(C|,^(e, T)/C'^g{e, T))/ In(Af) with respect to spreading gain M, when the path 
loss function is g{x) — ||a;||^" and e = 0.01. From the figure, we observe a similar Af^^^/a gain, even 
in the case of clustered transmitters. 



E=0.01,T=1,R.0. 25, 0.0.25,1 .1 




60 

M 



Fig. 9. \n{C*p„{e,T) /C}jg(e,T)) / \n{M) versus M for e = 0.01, Ap = 1 



V. Conclusions 

Previous work characterizing interference, outage, and transmission capacity in large random networks 
exclusively focused on the homogeneous Poisson point process as the underlying node distribution. In this 
paper, we extend these results to clustered processes. Clustering may be geographical, i.e., given by the 
spatial distribution of the nodes, or it may be induced logically by the MAC scheme. We use tools from 
stochastic geometry and Palm probabilities to obtain the conditional Laplace transform of the interference. 
Upper and lower bounds are obtained for the CCDF of the interference, for any stationary distribution of 
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nodes and fading. We have shown that the distribution of interference depends heavily on the path-loss 
model considered. In particular, the existence of a singularity in the model greatly affects the results. This 
conditional Laplace transform is then used to obtain the probability of success in a clustered network with 
Rayleigh fading. We show clustering the transmitters is always beneficial for large link distances, while 
the clustering gain at smaller link distances depends on the path-loss model. The transmission capacity 
of clustered networks is equal to the one for homogeneous networks. However, care must be taken when 
defining this capacity since clustered processes have two parameters to optimize over. We also show that 
the transmission capacity of clustered network is equal to the Poisson distribution of nodes. We anticipate 
that the analytical techniques used in this work will be useful for other problems as well. In particular the 
conditional generating functional are likely to find wide applicability. 
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Appendix 

A . Proof of Lemma 0" 

Proof: We first evaluate the asymptotic behavior of Q (^Fh ( g( -^) )) ■ ^'^^ ^(^^ ~ i^ gix-z) ) ■ have 
Mil v{x + u)f{u)dv\ = expf-c/ [l-Fh[^^ -]]f{u)du 



g{x + u — z) 

y 



- 1-c -) f{u)du (72) 

where (a) follows from the fact that exp(— x) = 1 — x + 0{x^) for x close to and (1 — Fh) for large 
y. By a similar expansion of exp, (j72p and the dominated convergence theorem, we have 



exp ~Xp / 1 — M{ / v{x + u)f{u)du)dx 



1- A. 



pC [ [ [I^fJ^-^ -)]/(w)d^.dx 



l-y-^^"XpC / [1-F;,(||x|r)]da; (73) 
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By change of variables and using limj,_oo[l - Fh{y)]y''''^ = [27, p.l98], we have 

/ [1 - Fh{\\xr)]dx = n r iy^/"dFh{u) 



We similarly have 



[ m( [ v{x + u)f{u)du] f{x)dx 



(a) 



Fh 



/(w)dw fix)dx 
g[x + u— z) J J 

-c [ f [1-Fh {y\\x + u- z\n] f{u)f{x)dudx 

-c[ [l-Fhiy\\x\n]{f*f){x + z)dx 

- cy-y- j [l-Fn (Ikir)] (/ * /) {-^ + z) dx 

-c{f*f){z)y-y'^ I [1-F,(||x|r)]dx 



(74) 



where (o) follows from the Lebesgue dominated convergence theorem (/ * / is a very nice function since / 
is a PDF). So we have Q (^Fh {^ g(JLz) ^ ~ B\y~'^l"-. For a Neyman-Scott cluster process, the second order 
product density is given by [16, p. 158], 

p(2)(,) = A^ + V(!lgp„„(„_i) 

TTC '-^ 

n=1 

where p„ is the distribution of the number of points in the representative cluster. /i(r)/7r denotes the density 
of the distribution function for the distance between two independent random points which were scattered 
using the distribution /(x) of the representative cluster. When the number of points inside each cluster 
is Poisson distributed with mean c, we have Y^=i'Pnf^kP' ~ 1) = c^- We also have /i(a;)/7r = (/ * /)(a;) 
Estimating ^p{%j) we have 

fiv) = — r / 9{x- z)p^'^>{x) / vdFh{v)dx 



\ r rVW^W 
= - / Ikll"" / i^dFh{v)d 

y Jm.^ Jo 



By change of variables, we have 



g r ry/g(x-z) 
+ — / (j(:r-z)ii(:r) / ;/dF,, (//)d.T 

T2 



Oi-2 Jo 



(75) 



(76) 
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For the term T2, 



Jo 



7ry2/a 

-2/Qr 



\\x\\''>v^^ ( + z j dx 

1 1 a: 1 1 > V 



= «i>5r.^'"di..(., ,77) 

So we have (/?(?/) ^ 02j/~^^"- Hence from Theorem [U we have F\{y) ^ 9iy~'^/"a,nd Ff{y) ~ {9i + ^2)2/^^^"- 



i?. Outage probability, in Poisson cluster process when the number of cluster points are fixed. 

In this subsection we derive the conditional Laplace transform in a Poisson cluster process, when the 
number of points in each cluster are fixed to be c S N and c > 0. We also assume that each point is 
independently distributed with density f{x). In this case the moment generating function of the number of 
points in the representative cluster is given by 

M{z) = z"" 

Using the same notation as in Section Hl-Bl and from (fTTj) and ([T3|) . we have 

= G{v) [ n v{x) [ nlid^y)fiy)dy 



G{v) f [ l[vix)^y{d4>y)fiy)dy 
~/r2 Jn ~- , 



'^^ G{v) [ ([ v{x - y)f ix)y ' f{y)dy (78) 

where (a) follows from the fact that the points are independently distributed and we are not counting the 
point at the origin. In this case G{v) is given by 



G{v) = exp |-Ap J^^ 1 - ^J^^ v{x + y)f{y)dy ) dx 
Hence the success probability (Rayleigh fading) is given by 



(success) = expl-Ap/" 1 - /3(i?, j/)%l /" j3{R,yf-^ f{y)dy (79) 



where 

iR2 1 + S^^g{x-y - z) 
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C. Outage probability of Nakagami-m fading 

Here, we derive the success probability when the fading distribution is Nakagami-m distributed. We also 
assume to G N and W = Q. The PDF of the power fading coefficient y ~ is given by 

From p5|l. we have 

P(success) = P ( T—-!^4^^^^—r-r > T] 



(success) = / — -(-) > y)dy (80) 



9{z)J Jo r(TO) \ nJ 



Using integration by parts we get, 



''(success) 



r( 



(a) (to - 1)! 

r(TO) 



1 pOC rri 

TO)/ r(TO,^,).P-(/,(.)<.) 



E M / e-^yy^dF'-\l^[z) < y) 



k=0 

(6) (^l)'" rf'^ 

= Yl -^^-^^l4.(z)i^)\s=Trn/ng(z) (81) 
fc=0 

where (a) follows from the series expansion of incomplete Gamma function when to is an integer and (b) 
follows from the properties of the Laplace transform and r(TO) = (to — 1)! when to is an integer. We also 
have 

1 



^hAsgix- z)) = 



Hence from Lemma [5J we have 



exp 



Xp l-ex.p{-cf3{s,z,y))dy 



eTq){-cP{s,z,y))f{y)dy 



(82) 



where 



/3{s,z,y) = 1 



-f{x)dx 



+ ^sg{x - y - z)y 

For integer to > 1, P(success) can be evaluated from (|8T|) and ([82]) . For to = 1, the probability evaluated 
from ([8T|l and ([82|l matches that of Lemma IH 
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D. Proof of Lemma^ 
Proof: 

Pp{Xpc) 



exp 
exp 



Xpc [ (3iR,y)dy + Xp [ (1 - exp[-c/3(i?, 2/)])d2; 
Xp / {c^(i?, y) - 1 + exp [ - c/3(i?, y)] }dy 



(83) 



Since 1 — exp(— ax) < ax, we have that v{R,y) > 0. We also have from the dominated convergence theorem 
and HI) 

^lim Ki?, 2/) - YTY^ - 1 + (-TTY^) > 

which is a constant. So using Fatou's lemma [27] (liminf J fn > J" (lini inf /„ ) . /„ > 0), we have 

hm ^^^^^'^^ 



(84) 



p"" Pfj^?^\\ " exp[-Ap lim / iyiR,y)dy] 
< exp[-Ap / lim v{R,y)dy] 
= exp[— ApOo] — 



E. Proof of Lemma\S[ 

Proof: From ([55l) . the probability of success is 



(success) = Pp(Apc)exp ApC / l3{R,y)'q{c,R,y)dy P2{R,c) 

V ' T2 

Ti 



(85) 



where ri{c,R,y) — v{c(3{R,y)) and v{x) = (exp(— .t) — 1 + x)/a; an increasing function of x. From Young's 
inequality [27, Sec. 8.7] we have l3{R,y) < min{l, sup{/(x)}i?2T2/"C(a)}. Hence 

r]{c,R,y) < z^(cmin{l,sup{/(x)}i?2T2/"C(a)}) 
With a slight abuse of notation, let 77(0, R) = i/(cmin{sup{/(a;)}i?^r^/"C(a), 1}). Hence 

Ti < exp[ApC / P{R,y)rjic,R)dy] 
= exp[ApCr2/"i?2r;(c, R)C{a)] 
Also observe that 7j{c, R) ^ R'^. So Ti < exp(i?4). 



~ 00 
T2 = l-cf3{R,y) + cp{R,y)J2' 



("1)" 



c/3(i?,y))"-V(y)dy 



< 



[1 - c/3(i?, y) + c/?(i?, y)r,{c, R)]fiy)dy 
l-[l-77(c,i?)]c / /3(i?,y)/(y)dy 



(86) 



(87) 
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If one considers x and y as identical and independent random variables with density functions /, we then have 

/g2 P{R,y)f{y)dy — E[ — ^(-BTTr 17^] ■ < k < 1 be some constant. Using the Chebyshev inequality 

we get 



E 



1 



l + &\x-y-R\\ 



> kP 
= kP 



1 



I + ^\\x - y - RW^ 
1 



> K 



y-R\\<i--iy/°'RT^^°' 



(88) 



The PDF of z = a; — J/ is given by (/ * f){z), since y is rotation- invariant. Choosing k = T/(l + T) we have 

T 



1 + T 



{f*f){x)Ax 



> 



T 



1 + T 



B(R,R) 



B(R,R) 



\x\\dx 



1 + ^ Jb{i,i) 



ja;||dx 



(89) 



C2 



So we have 



P2 < 1 - [1 - ?7(C, i?)]i?'C2 

< l-R^ + R^ 



(90) 



Also we have Ti < exp(i?4) < i + l.01i?4, go we have P2T1 <1-R^ + R^ - l.Oli?^ + l.Oli?^ < 1 for smaU 
R 0. Hence for small R we have ¥{success) < Pp(Xpc). 



F. Proof of Lemma \ll\ 

Proof: We find C„(e, T) and hence upper bound the transmission capacity. We have from the derivation 
of Lemma [5] 



P(A, T) < PpiXpc) exp[XpCf3irj{c, R)]P2{R, c) = P„(cAp, T) 



(91) 



where r]{c,R) — (exp(— c/3) — 1 + cj3)/c(}. With — {\pC, P^{\pC,T) > 1 — e}, it is sufficient to prove 
sup^u < Cp{e,T). Also observe that Pu{cXp,T) — > as c — > 00 independent of Ap. Hence we can assume c 
is finite for the proof. We proceed by contradiction. 

Let sup v4„ > Cp{e,T). Hence there exists a S > 0, Xp > 0, c> such that XpC = (Cp(e, T)/ {l—e))+S S A^- 
At this value of XpC we have 

l-e<Pu(cAp,T) = il-e)Pp{5)exp[r^{c,R){Hj^J+Sl3i}]P2iR,c) 

= {l-ey-^^'^''^Pp{S{l^rj{c,R)))P2{R,c) (92) 
« ■' 

Ti 
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From the derivation of Lemma [SJ we have P2{R,c) < 1 — [1 — iri{c, R)\ck, with equahty only when c = 0. 
Hence we have 

Pu{c\p,T) < ri(l-e)i-"(^"'^)(l- [l-rKc,i?)]cA:) (93) 
Since ''^p(^^^^(^^^) < -^^e have r]{c, R) < c$/{l + c$). Using the upper bound for 77(0, R) 

0(3 



Pu{c\p,T) < Ti(l-e)(l-e) 1+^/5 1 _ 



1 - 



l + c(3 



cPp{T) 



= T,,i-,)(i-.rrf&(i-fS) 

V 1 + c/3 / 



(94) 



T2 



Using the inequahty I — ay < (1 — b)^ , b < 1 ~ e ",2/>0, substituting y — b = e,a ^ PiT), we get 

T2 < 1. Hence we have 



PuicXp,T) < (l-e)Pp(,5(l-77(c,i?))) 



(95) 



So if (5 > 0, and c finite, we also have Pp{S{l — 77(5, R))) < 1. So we have a contradiction from ([5^ and ([55)1 . 
Hence there exists no such S and hence supA„ < Cp{e,T). We can achieve Cu{e,T) — Cp{e,T), by using 
Xp = n^i^j^. — 1/71 for rt very large. As n 00, Pu{cXp,T) Pp{cXp,T) . ■ 
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